Intrinsic uncertainty on the nature of dark energy 
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We argue that there is an intrinsic noise on measurements of the equation of state parameter 
W = p/ p from large-scale structure around us. The presence of the large scale structure leads to an 
ambiguity in the definition of the background universe and thus there is a maximal precision with 
which we can determine the equation of state of dark energy. To study the uncertainty due to local 
structure, we model density perturbations stemming from a standard inflationary power spectrum 
by means of the exact Lemaitre-Tolman-Bondi solution of Einstein's equation, and show that the 
usual distribution of matter inhomogeneities in a ACDM cosmology causes a variation of w - as 
inferred from distance measures - of several percent. 
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Introduction A key quantity to characterize the na- 
ture of the dark energy is its equation of state parameter 
w = p/ p. Current and future cosmological observations 
try to measure w ever more accurately, and the power of 
dark-energy missions is judged by the minimal error that 
they can achieve on the dark-energy equation of state 
parameter w = p/p. This is for example the basis of 
the Dark Energy Task Force (DETF) [1] Figure of Merit 
(FoM), which is given by the determinant of the Fisher 
matrix for the parameters Wq and w a in a linear param- 
eterization of the equation state, w(a) = wq + (1 — a)w a . 
An important question in the context of dark-energy re- 
search is then whether there is a natural limit for the 
precision with which w can be measured, or whether one 
can in principle determine w(a) to an arbitrary precision. 

In this letter we argue that the matter fluctuations that 
are always present in the universe provide such a limit, 
and we determine the unavoidable variation of w(a) as 
expected in the ACDM concordance model (for which in 
principle w = —1). Those variations appear because we 
always measure any observable quantity in the true per- 
turbed universe, even if we consider "background" quan- 
tities like the luminosity distance (see e.g. [2] and ref- 
erences therein), and they remain significant even when 
averaging over angles [3] . This is a direct manifestation of 
the "fitting problem" [4], i.e. the attempt to fit a homoge- 
nous and isotropic FLRW model to a lumpy universe [5- 
7] rather than directly modeling the inhomogeneities [8]. 
If, on one hand, fluctuations in e.g. the luminosity dis- 
tance allow us in principle to obtain additional cosmo- 
logical information (see e.g. [9-13]), on the other hand 
they result in an intrinsic noise in the determination of 
cosmological parameters. Indeed, although the pertur- 
bations in the metric are small, only about 10 -5 , they 
can be amplified when going to quantities that involve 
derivatives like w(a), as demonstrated in e.g. [14]. 



The outline of this letter is as follows: first we define 
the probability of inhomogeneities. Then we explain how 
we model said inhomogeneities. Next we show the noise 
caused by this structure on the evolution of the dark- 
energy equation of state inferred by an observer who ig- 
nores the inhomogeneity. Finally we obtain the error 
in the cosmological parameters measured by an observer 
who fits the luminosity distance under various assump- 
tions. 

Probability of a local structure The root mean 
square of the density perturbation in a sphere of radius 
L around any point today in a Gaussian density field 
is [15] 
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where P(k) is the matter power spectrum today as a 
function of wavenumber fc, and ji is the spherical Bessel 
function of the first kind. 

Given some inhomogeneity with a mass M(L) inside 
a radius L, one can define the average density perturba- 
tion Sq = M(L)/M(L) — 1 relative to the homogeneous 
background which predicts a mass M(L) inside the same 
radius L. Then the probability of having such a structure 
is [16], 
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If for simplicity we model the inhomogeneity spherically 
symmetric, that is, we presume that the central observer 
averages observations over directions, the mass is given 
by M(r) = 4-7T J Q dr^/—gp rn {r) with g the determinant of 
the metric. 
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Model for local inhomogeneity To model the in- 
homogeneity averaged over angular directions, we adopt 
the spherically symmetric Lemaitre-Tolman-Bondi solu- 
tion [17-19] including a cosmological constant A (ALTB, 
see e.g. [20-23]), for which we can compute all distance 
measures exactly. 

The ALTB metric in the comoving and synchronous 
gauge can be written as (using units for which c = 1) 



ds 2 = -dt 2 



ap,r) 
1 — k(r)r' 



-dr 2 + al(t,r)r 2 dn 2 , (3) 



where the longitudinal (ay) and perpendicular (a_i_) scale 
factors are related by ay = (a±r)', and a prime denotes 
partial derivation with respect to the coordinate radius 
r. In the limit k —> const., and a± = au we recover the 
FLRW metric, but in a LTB metric the curvature k(r) is 
a free function and in general is not constant. 

The two scale factors define two different Hubble rates: 



0_L 

a± 



H\\(t,r) 



(4) 



The analogue of the Friedmann equation in this space- 
time can be written in a familiar form, 
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where we adopted the gauge fixing a± — 1. However, 
the density parameters are now also functions of r, 
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so as to satisfy f2 m (r) + f2&(r) + ^a(^) = 1- The 
free function m(r) is related to the local matter den- 
sity 8ttG p m (t,r) — (mr 3 )' / 'a\\a\r 2 . The matter den- 
sity obeys the usual conservation equation p' m + (2H± + 

H \\)Pm = 0. 

Finally, time t and radius r as a function of redshift 
z are determined on the past light cone of the cen- 
tral observer by the differential equations for radial null 
geodesies, 
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with the initial conditions i(0) = to an d ^(0) = 0. The 
area {(1a) and luminosity (dr,) distances are given by 

d A (z) =a x (t(z),r(z)) r(z), d L = (l + z) 2 d A . (8) 

Density profile The age of the universe is a function 
of (t, r) and is obtained by integrating the Friedmann 
equation (5) from the big-bang time tbb(^) to time t: 
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Eq. (9) relates the three free functions <bb, k(r) and m(r), 
so that density of the dust field in the ALTB model is 
specified by two free functional degrees of freedom, where 
we choose k(r) and t^(r). Any radial dependence of 
tbb( r ) is directly related to a decaying mode in the mat- 
ter density field [24, 25]. By choosing tbb( r ) = decaying 
modes are absent, in agreement with the standard sce- 
nario of inflation. 

We parameterize the curvature function with the 
monotonic profile 



k(r)=k b + (k c -k b ) P 3 (r/r b ) 
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where r b is the comoving radius of the spherical inhomo- 
geneity and P3 is the function 



Pn{x) 
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for n = 3. The function P n (x) interpolates from 1 to 
when x varies from to 1 while remaining n times 
differcntiablc, which implies that that k(r) is C" every- 
where. We choose n = 3, such that the metric is C 2 
and the Riemann curvature is C . For r > r b the cur- 
vature profile equals the curvature k b of the background 
such that there the metric reduces exactly to the ACDM 
model. The central under- or over-density, determined by 
the curvature k c at the center, is automatically compen- 
sated by a surrounding over- or under-dense shell. We 
adopt the conservative approach of using a compensated 
density profile so as not to alter the background metric 
of the universe, which otherwise would be FLRW only 
asymptotically. The radius L of the inhomogeneity that 
is used in Eq. (2) is the radius at which the central over- 
or under-density has the transition to the surrounding 
compensating under- or over-dense shell. The radius L 
is hence smaller than the radius r b which defines the ra- 
dius of the total LTB patch, including both the central 
perturbation and its compensating shell. 

In summary, the local structure is parametrized by the 
radius of the boundary r b and the central curvature k c . 
For any choice of these parameters, and for a ACDM cos- 
mology (given by the background matter density, the cur- 
vature parameter k bl the present-day background Hubble 
rate and a specific P(k)), we can compute the probabil- 
ity of the existence of such a structure using Eq. (2) and 
compute the luminosity distance from an object to the 
central observer as a function of redshift. 

FLRW Observer's w(z) Following [14] one can - 
given a luminosity distance-redshift relation in a homo- 
geneous universe described by the FLRW metric - com- 
pute what the underlying w(z) of the dark-energy fluid 
is. There is an exact relation between w(z) and the first 
and second derivatives of the luminosity distance with re- 
spect to redshift and two more parameters, and f2 m . 
Therefore, the observer can derive w(z) if he/she knows 
the latter two parameters from other observations, and 
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Figure 1. Matrix representation of the contribution of matter 
perturbations with a radius 0.1 (i — 1) < Zb < O.li (the index 
i = 1..12 labels rows) to the dispersion a w r x ) at Zj = 0.05 + 
O.lj (the index j = 1..12 labels columns). 



deduces the derivatives of the luminosity distance from 
SN observations. In the scenario studied here, w = — 1. 
However, the inhomogeneities come into play and modify 
the luminosity distance-redshift relation. Consequently, 
an observer that erroneously assumes that the metric sur- 
rounding him/her is FLRW will in fact see a redshift- 
dependent w [8, 26-28]. 

In order to study this fundamental variation of the 
reconstructed w(z), we first set the fiducial flat ACDM 
model to the WMAP7+LRG best-fit cosmology [29]. We 
then sample the {rb,k c } parameter space by building a 
Markov chain using Eq. (2) as likelihood, but restricted 
to a certain range in redshifts 0.1 (i — 1) < Zb < O.li 
for i = 1..12, where Zb = z(rb) is the apparent redshift 
at which an observer sees the radius rb (see Eq. (7)). 
Next we compute from Eq. (8) the corresponding lumi- 
nosity distances in the various realizations. We finally 
derive using Eq. (3) of [14] the (apparent) w(z) that an 
FLRW observer would infer. In this procedure we let the 
FLRW observer fix il m to the fiducial value as it can- 
not be determined from cosmological observations alone 
when allowing for an arbitrary equation of state of the 
dark energy [30]. The curvature could in principle be 
constrained by combining measurements of the distances 
with measurements of the expansion rate H(z), but for 
simplicity we also fix fi^ to the fiducial (zero) value. The 
result of this analysis gives the variance in w{z) at the 
nodes Zj = 0.05 + O.lj for j = 1..12 induced by structures 
falling in the redshift bin 0.1 (i — 1) < Zb < O.li, which is 
shown in Fig. 1. In Fig. 2 we show examples of typical 
uu(z) evolutions as seen by the FLRW observer, affected 
each by one structure of an arbitrary size. 

Fitting the luminosity distance However, usually 
one does not derive a fully general w{z) but fits a parame- 
terized model to the distance data. Typical examples are 
a constant w or the linear model used for the DETF FoM 
mentioned in the introduction, w(a) = Wq + (1 — a)w a . 
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Figure 2. Examples of apparent evolution of the dark-energy 
equation of state that an FLRW observer would deduce from 
observations in a flat ACDM universe endowed with a local 
inhomogeneity coming from a post-inflationary gaussian den- 
sity field. 



As the impact of the local structure is strongest at low 
redshift, the variance of the fitted parameters will de- 
pend on the way the weight of distance data depends 
on redshift. For the purpose of our letter we choose a 
specific redshift distribution modeled to resemble the ex- 
pected supernova distribution of the Dark Energy Survey 
(DES). Specifically, we weight the 12 redshift bins in the 
range < z < 1.2 by using the binned rms scatter Ubin 
of the Simulated DES Hybrid 10-field Survey reported in 
the third column of Table 14 of [31] . 

The basic approach is then as above: We fix again a 
fiducial cosmology with the same parameters as in the 
previous section (WMAP7+LRG) and sample again the 
parameter space describing the inhomogeneities in sepa- 
rate redshift bins 0.1(i — 1) < Zb < O.li using Eq. (2) as 
likelihood. We then fit the parameterized model to the 
data and determine the best-fit parameters 9* by mini- 
mizing the following x 2 '- 



x 2 m 
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where Zj = 0.05 



O.lj, mhom and mi n h are the distance 
moduli of the homogeneous ACDM model and of the in- 
homogeneous model, respectively (with the latter playing 
the role of the data in our context), and we marginalize 
analytically the likelihood oc e~ x I" 1 over an unknown off- 
set. This time we can optionally also vary and fi m as 
the parameterized model breaks the degeneracies. 

In Fig. 3 we show the dispersion induced by a sin- 
gle structure of size Zb < 0.5 on the fit parameters 
{£l m ,Wo,w a } when the observer assumes that Q], = 0. 
Since the effect does not change much for different struc- 
tures up to Zb = 0.5, it is safe to consider the average 
variance induced by such structures. We compute the 
expected variance for such structures by combining the 
MCMC chains for the size bins up to Zb < 0.5 by thin- 
ning them down to having the exact same number of 
points and then merging the resulting thinned chains. 



4 



o 0.1 



1 1 


1 1 

o 


w ° FEE! 




o 
























■ 

























0.1 0.2 0.3 0.4 0.5 

z b 

Figure 3. Dispersion induced by a single structure of size 
Zb < 0.5 on the fit parameters {fi m , wo, ifa} when Qk = is 
assumed. The last column in Table I gives the corresponding 
values when one averages the effect of these structures. 

This leads to a single chain in which each size bin con- 
tributes with equal weight. We list in Table I the nu- 
merical values of the dispersions for the four different 
parameterized models. The variances listed in Table I 
show that the density perturbation around us, stemming 
from a standard inflationary spectrum and of unknown 
density, adds an extra uncertainty to these parameters. 
Unless a measurement of the local perturbation becomes 
possible in the future, we can never measure these pa- 
rameters at a higher accuracy then the listed variances. 





{O ro ,Ofe} {fl m ,W } {Q. m ,Qk,W } {fl 


■n, Wo, W a } 




0.013 0.0043 0.019 


0.018 




0.029 0.067 






0.020 0.057 


0.025 






0.18 



Table I. Intrinsic la uncertainties on fitted parameters for 
four different FLRW models. In each model all remaining cos- 
mological parameters are fixed to the WMAP7+LRG best-fit 
cosmology. The free parameteres are fitted to the inhomoge- 
neous distance by minimizing the \ 2 °f Eq. (11). The listed 
dispersions are the standard deviation of the posterior distri- 
bution on the given parameter, marginalized over the other 
free parameters. 



then A min = TiVdetCUnAx 2 = Tr<T Wo a Wa y/l - p 2 Ax 2 
where A\ 2 — 5.99 for a 95% c.l. contour. We find 
l/x/detC min = 581 and FoM max = 31. 

Conclusions We have estimated in a conservative 
way the intrinsic uncertainty in the reconstruction of the 
dark-energy equation of state by means of distance mea- 
surements. Since we observe only one universe, this un- 
certainty will show up as a bias in distance measure- 
ments, and we propose that the scientific community 
use the results of this letter so as to include this ex- 
tra source of error in their analysis. In particular, we 
give the covariance matrix C m ; n for the linear model 
w(a) = w + (1 - a)w a used for the DETF FoM, which 
can be easily convolved with any other posterior distri- 
bution constraining the parameters wq and w a . Since the 
large-scale structure limits strongly the power of distance 
measurements as a probe of the nature of the dark en- 
ergy and of the curvature of the universe, one may need 
to use data e.g. from galaxy surveys, weak lensing mea- 
surements or from the integrated Sachs- Wolfe effect to 
reduce its impact at least partially. 

This analysis can be extended in a number of ways. For 
example, one may model the local inhomogeneity using 
non-sphcrically symmetric and non-compensated density 
profiles. Even more important would be the inclusion 
of non-local structures like superclusters and filaments 
which may increase the noise through lensing. According 
to results from second-order perturbation theory [28, 32] 
there is a scatter of about <r M w 0.02 or more in the 
luminosity distance even for z > 0.4 which can lead to a 
percent-level variation for w also at high redshifts. 
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Regarding the FLRW model where {£! m , wq, w a } were 
left free, we observe that - based on these results - an 
experiment like DES can never determine wq to a pre- 
cision better than 2.5%, and w a better than 18%. The 
Figure of Merit is defined as FoM= 1/A where A is the 
area bounded by the 95% c.l. contour on the {u>o,u> a } 
plane. The covariance matrix for {wq, w a } (which is the 
inverse of the Fisher matrix) for this intrinsic noise on 
the dark-energy equation of state is: 

C mi n=( P(Tw f w A, (12) 

\pcr Wa a Wa a Wa ) 

where a Wo and o~ Wa are given in the last column of Ta- 
ble I and the correlation is p = —0.924. The area is 



[1] A. Albrecht, G. Bernstein, R. Cahn, W. L. Freedman, 
J. Hewitt, et ai, (2006), arXiv:astro-ph/0609591 [astro- 
ph]. 

[2] C. Bonvin, R. Durrer, and M. A. Gasparini, Phys.Rev. 

D73, 023523 (2006), arXiv:astro-ph/0511183 [astro-ph]. 
[3] N. Mustapha, B. Bassett, C. Hellaby, and G. Ellis, 

Class. Quant. Grav. 15, 2363 (1998), arXiv:gr-qc/9708043 

[gr-qc]. 

[4] G. Ellis and W. Stoeger, Class.Quant.Grav. 4, 1697 
(1987). 

[5] E. W. Kolb, V. Marra, and S. Matarrese, Gen.Rel.Grav. 

42, 1399 (2010), arXiv:0901.4566 [astro-ph. CO]. 
[6] T. Clifton and J. Zuntz, Mon.Not.Roy.Astron.Soc. 400, 

2185 (2009), arXiv:0902.0726 [astro-ph.CO]. 
[7] P. Bull and T. Clifton, Phys.Rev. D85, 103512 (2012), 

arXiv:1203.4479 [astro-ph.CO]. 



5 



[8] V. Marra, M. Paakkonen, and W. Valken- 

burg, Mon.Not.Roy.Astron.Soc. in press (2013), 

arXiv:1203.2180 [astro-ph.CO]. 
[9] P. Valageas, Astron.Astrophys. 354, 767 (2000), 

arXiv:astro-ph/9904300 [astro-ph] . 
[10] S. Dodelson and A. Vallinotto, Phys.Rev. D74, 063515 

(2006), arXiv:astro-ph/0511086 [astro-ph]. 
[11] D. Sarkar, A. Amblard, D. E. Holz, and A. Cooray, 

Astrophys.J. 678, 1 (2008), arXiv:0710.4143 [astro-ph]. 
[12] C. Bonvin, R. Durrer, and M. Kunz, Phys.Rev.Lett. 96, 

191302 (2006), arXiv:astro-ph/0603240 [astro-ph]. 
[13] L. Amendola, K. Kainulainen, V. Marra, and M. Quar- 

tin, Phys.Rev.Lett. 105, 121302 (2010), arXiv: 1002. 1232 

[astro-ph.CO]. 

[14] C. Clarkson, M. Cortes, and B. A. Bassett, JCAP 0708, 
011 (2007), arXiv:astro-ph/0702670 [ASTRO-PH]. 

[15] E. W. Kolb and M. S. Turner, Front. Phys. 69, 1 (1990). 

[16] W. Valkenburg, V. Marra, and C. Clarkson, (2012), 
arXiv:1209.4078 [astro-ph.CO]. 

[17] G. Lemaitre, Ann.Soc.Sci.Bruxelles 53, 51 (1933). 

[18] R. C. Tolman, Proc.Nat.Acad.Sci. 20, 169 (1934). 

[19] H. Bondi, Mon.Not.Roy.Astron.Soc. 107, 410 (1947). 

[20] B. Sinclair, T. M. Davis, and T. Haugbolle, Astrophys.J. 
718, 1445 (2010), arXiv: 1006.0911 [astro-ph.CO]. 

[21] V. Marra and M. Paakkonen, JCAP 1012, 021 (2010), 



arXiv:1009.4193 [astro-ph.CO]. 
[22] W. Valkenburg, Gen.Rel.Grav. 44, 2449 (2012), 

arXiv: 1104. 1082 [gr-qc]. 
[23] A. E. Romano and P. Chen, JCAP 1110, 016 (2011), 

arXiv: 1104.0730 [astro-ph.CO]. 
[24] J. Silk, A&A 59, 53 (1977). 

[25] J. P. Zibin, Phys. Rev. D78, 043504 (2008), 
arXiv:0804.1787 [astro-ph]. 

[26] A. Shafieloo, V. Sahni, and A. A. Starobinsky, Phys.Rev. 
D80, 101301 (2009), arXiv:0903.5141 [astro-ph.CO]. 

[27] G.-B. Zhao, R. G. Crittenden, L. Pogosian, and 
X. Zhang, Phys.Rev.Lett. 109, 171301 (2012), 
arXiv: 1207. 3804 [astro-ph.CO]. 

[28] I. Ben-Dayan, M. Gasperini, G. Marozzi, F. Nugier, 
and G. Veneziano, Phys.Rev.Lett. 110, 021301 (2013), 
arXiv:1207.1286 [astro-ph.CO]. 

[29] E. Komatsu et al. (WMAP Collaboration), Astro- 
phys.J. Suppl. 192, 18 (2011), arXiv: 1001.4538 [astro- 
ph.CO]. 

[30] M. Kunz, Phys.Rev. D80, 123001 (2009), arXiv:astro- 

ph/0702615 [astro-ph]. 
[31] J. Bernstein, R. Kessler, S. Kuhlmann, R. Biswas, 

E. Kovacs, et al, Astrophys.J. 753, 152 (2012), 

arXiv:1111.1969 [astro-ph.CO]. 
[32] I. Ben-Dayan, M. Gasperini, G. Marozzi, F. Nugier, and 

G. Veneziano, (2013), arXiv: 1302.0740 [astro-ph.CO]. 



